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We continue the program started in |7| and explain the statistics of the excitations for the 
generalizations of the paired states in the quantum Hall effect in terms of the parafermion statistics. 
We show that these excitations behave as combinations of bosons and parafermions. That generalizes 
the prior treatment of the paired (PfafRan) state where the excitations behave as combinations of 
bosons and fermions. We explain what it means, from a quantum mechanical point of view, for a 
particle to be a 'parafermion' rather than a boson or a fermion and work through several explicit 
examples. The resulting multiplets coincide exactly with the angular momentum multiplets found 
numerically for the k + 1 particle interaction Hamiltonian on a sphere. We also present a proof that 
the wave functions found in are indeed the correlation functions of the parafermion conformal 
field theory. 



I. INTRODUCTION 

The quantum Hall plateaus on the first excited Landau levels appear to be different from those of the lowest Landau 
level. The most striking difference is the 5/2 plateau |l| for which there is no analog in the lowest Landau level. 

To explain the 5/2 plateau, several trial wave functions were proposed. One of them, the Haldane-Rezayi state [g), 
assumes that the electrons in it are spin unpolarized. Yet another trial wave function, the Pfaffian assumes that 
electrons are polarized, just as for the ordinary Laughlin trial wave functions Recent numerical evidence in 
the absence of the Landau level mixing, supports the Pfafnan as a trial wave function for the 5/2 plateau. 

In the paper by N. Read and one of us a further generalization of the PfafRan state was proposed. Moreover, 
it was conjectured, and supported by the numerical evidence, that these generalizations may be better candidates 
for other plateaus on the first excited Landau levels than the Laughlin state and its hierarchical ||,|]f]] or composite 
fermion |n| generalizations. 

An unusual feature of the generalizations of the Pfaffian state given in was that the trial wave functions were 
found explicitly for nonsimple fractions. That was done with the help of conformal field theory. 

Conformal field theory discovered in jll[] is essentially a method to solve various scale invariant 1+1 dimensional 
quantum field theories exactly. It was since proved extremely useful for understanding various 1 dimensional quantum 
and 2 dimensional classical statistical mechanics systems. Its relevance for the quantum Hall systems was first discussed 
in • 

The relevant conformal field theory for the generalization of the Pfaffian is the so-called parafermion conformal 
field theory [ p~2| . The parafermions, the direct generalizations of the fermions, play an important role in describing 
the critical points of the Z^ invariant statistical mechanics systems. A particularly well known example of the Z 2 
invariant system is the Ising model. The Ising model can be described by a Majorana fermion. A Zk invariant system 
must be described, on the other hand, by Z/j-parafcrmions. 

A Majorana fermion ^p(z) can be used to generate the Pfaffian trial wave function, as was discussed in 

^Pfaffian = (lp(zi)lp(z 2 ) ■ ■ ■ 1p(z N )) j [ (Zi ~ Zj) 2 (1) 

i<j<N 

The filling factor of this trial wave function is 1/2. Moreover, the so-called spin fields of the fermion conformal field 
theory (the order parameter of the Ising model) a can be used to generate the excitations of such a system, 

*p X faffia d „(^1^2, ■ • .,V2L) = ((?{Vl), <*(»»)> • • ■ , * fei), VK^), . . . ,^{z N )) JJ(ffc - Zj)* JJfo - Zjf . (2) 

On the other hand, as was discussed in H], for each k there are k — 1 parafermions, denoted ipi, with dimensions 
hk = l(k — l)/k. In it was proposed to use the first of them to generate a parafermion wave function, 
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*jL,=^(^iw-^w n (*i-*o) M+i (3) 

i<j<AT 

where M has to be taken odd or even integer depending on whether the particles are fermions or bosons. The filling 
factor of this wave function is given by v — k/(kM + 2). 

It was further shown in |7j that this wave function should be an exact ground state of a system of N bosons in a 
magnetic field in the presence of the fc-body interaction Hamiltonian 

H = V J2 S 2 ( Zil -z i2 )5 2 (z i2 -z i3 )...S 2 (z ik -z ik+1 ), (4) 

il<i2<---<ik + l 

thereby generalizing the Pfaffian case for which the interaction Hamiltonian was (Q) at k = 2. We do not need to 
consider fermions because their wave function can be obtained from that of bosons by a simple Jastrow factor. 

To find the wave function directly from (^) is a difficult task. However, it was shown in Q that the following explicit 
construction has the correct properties (it is symmetric and vanishes whenever k + 1 particles coincide), which we 
reproduce for future reference. 

To write down ^para we need to break the coordinates of N electrons into clusters of k (N should be divisible by 
k). For each pair of distinct clusters, say z±, . . . , z k and Zk+i, . . . , z 2k , we define factors \ by 

Xi, 2(2:1, ■ ■ ■ ,Zk;Zk+i, ■ ■ - ,Z2k) = (zi ~ z k+1 )(zi - z k+2 ){z 2 - z k+2 ){z 2 - z k+3 ) ...{z k - z 2k )(z 2k - z k+ i). (5) 

The subscript 1,2 of xi,2 labels the clusters (the first, starting with Z\ and the second, starting with z k ). 
The wave function ^pala is defined in terms of these as 

*pala = II Xr+l,s+l(2 : P(fcr+l),--- I 2 : p( fe ( r+1 ));Zp( fe;s+1 ),...,Zp( fc ( ;5+1 ))). (6) 

P 0<r<s<N/k 

It is not hard to see that for k = 2 this wave function reproduces the Pfaffian for bosons at v = 1. 

It was also observed in that the quasihole excitations above (||) can be described by the insertion of the spin 
fields a i of the parafermion conformal field theory. While some of the excitations for the Pfaffian state were found 
explicitly in [fl3||, no explicit wave functions for the excitations of the parafermion states have yet been found. 

The purpose of this paper is to continue the work begun in 0. In particular, we explain the numerically observed 
degeneracies of the excitations of the parafermion states. It turns out that to reproduce the numerically observed 
degeneracies we need to assume that these excitations obey the parafermion statistics. 

Here we would like to emphasize that although the parafermions have been shown in jl5) to obey the Haldane 
exclusion principle in the statistical sense which is a valid description for large quantities of parafermions, for 
small number of parafermions we need to know the concrete rules which govern the combinatorics of parafermions. 
These rules asymptotically approach the exclusion statistics rules as the number of parafermions becomes large. These 
rules are described below and by themselves do not have much to do with Haldane exclusion statistics. 

Our finding that the excitations of the parafermion states behave as a combination of bosons and parafermions is 
a direct generalization of the well known fact that the excitations of the Pfaffian state behave as a combination of 
bosons and fermions. 

We also show that the wave function (^) is indeed given by the parafermion correlation function (j^) , thus confirming 
the conjecture of M. 



II. EXCITATIONS OF THE Z 3 -PARAFERMION STATE ON A SPHERE 

In this section we are going to explain how to calculate the degeneracy of the excitations of the parafermion states. 
To do that, we will have to construct what could be called parafermion quantum mechanics and learn how to add 
angular momenta of the quantum mechanical particles obeying parafermion statistics. 

It was shown in M that if you put N particles on a sphere with a magnetic monopole in the center, turn on k + 1 
particle interaction^) and then adjust the total flux of the magnetic field through the surface of the sphere to be 
equal to = v^ 1 N — [M + 2) with v being the filling factor, v = kj (Mk + 2) and M being either odd (for fermions) 
or even (for bosons) nonnegative integer, you discover that the electrons settle into one unique zero energy ground 
state with the wave function 
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Kj 

with ^(^i, . . . , Zn) given by (^|). The total angular momentum of that state is equal to 0. 

Now if you start increasing the flux by units of one flux quanta, you discover that the zero energy eigenstates of 
([|) are degenerate. These states can all be grouped into the angular momentum multiplets. 

The reason for this is more or less clear. By increasing the flux by 1, we create quasiholes. There is more than one 
way of creating those quasiholes which explains the degeneracy of states at a higher flux. 

As was explained in Jt| the elementary quasiholes at arbitrary k are not Laughlin quasiholes. Rather, they are 
quasiholes which carry flux X/k and we create k of them at once. The wave functions with quasiholes can be found 
with the help of conformal field theory, by inserting the fields a\ of the parafermion conformal field theory Jl^,^[. 
However, we are not going to be interested in the explicit wave functions. 

What we want to explain in this section is the numerically observed degeneracy of the states with quasiholes on a 
sphere and how exactly one could generate all the degeneracies. We will show how the correct angular momentum 
multiplets can be obtained by putting parafcrmions into orbitals on a sphere and combining their angular momenta 
in a way consistent with their fractional statistics. That generalizes the prior treatment of the Pfaffian state in jlj]. 

Now we would like to present, first without explanation, the rules for finding the degeneracies of the excitations 
above the k = 3 state, that is, for Z3-parafermions. 

It was found in [Q, by matching the numerical data, that, as we increase the flux, the degeneracy of states we get for 
the Z% parafermions is given by the following formulas. The number of excitations at the excess flux 1 (3 quasiholes, 
each carrying I /3 of the flux quantum) is given by a binomial coefficient 



N/3 + 3 



(8) 



with N being the total number of the electrons. N should be divisible by three. This coefficient has a simple 
interpretation as the number of ways you can put 3 bosons into N/3 + 1 orbitals. The reason why the excitations at 
the excess flux I behave as bosons was explained in JtJ and is essentially the same for the parafermion states as for 
the Pfaffian state JTJ. 

Moreover, following pi , we can assign the angular momentum quantum numbers to these states. For this purpose 
we interpret the N/3 + I states as the orbitals on a sphere, in the multiplet of the angular momentum L = N/6. 
When we put the bosons on this sphere, we generate the angular momentum multiplets by combining their angular 
momenta while keeping their wave function totally symmetric. 

While we could refer to any standard quantum mechanics textbooks for the rules on how to combine the angular 
momenta of bosons, we can recast these rules into the following simple form. Let us visualize ./V/3 + I orbitals as 
boxes. Each box has a number L z assigned to it which varies from —N/6 to N/6 by steps of f. We are allowed to 
put as many bosons as we wish into each particular box. 

For example, for N = 6 we have 3 boxes with L z = —I, L z = 0, and L z = 1. Fig. [I] shows one possible way to put 
3 bosons into 3 boxes, by putting 2 of them into L z = 1 box and one into L z = — 1 box, thus giving the total L z = 1. 
Of course, there are (5) = 10 ways to put 3 bosons into 5 boxes. Out of these 10 ways, the angular momentum 
projections L z — ±3 or L z = ±2 can be obtained in 1 way each, while the angular momentum projections L z = ±1 or 
L z = can be obtained in 2 different ways each. We immediately recognize the sum of the L — 3 and L = I angular 
momentum multiplets. And indeed, this is confirmed by the numerical data. 



FIG. 1. One possible way to put 3 bosons into 3 boxes, realizing M for N = 6 



As the flux is increased we observe many more excited states. It was found in Q that the number of states at the 
excess flux 2 is given by 



N/3 + 6 
6 



N/3 + 5 
6 



fN/3 + 4 



(9) 
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and at the three excess flux quanta 



( w/3 9 +9 ) + i°H + >'»( w/3 9 +7 )' (io) 

The higher fluxes turned out to be harder to analyze as there is less numerical data available for them. 
We observe that the number of states at these flux quanta is reminiscent of the corresponding formula for the 
Pfaffian which was found in Hi to be 



/n\ f{N - F)/2 + 2n 

F,(-1) F = (-1)« ,F<N ^ ' ^ 



(11) 



with n being the number of the excess flux quanta. The extra (^,) was interpreted as the number of ways one can put 
F fermions into n boxes. By analogy, we can write down a formula which generalizes ( fTT|) to the case of the k = 3 
parafermions, 

£ {;}{ (N ~T +3n ) <«> 

F>0,exp(22fi) = l,_F <N 

where { ^ } is the number of ways you can put F i^-parafermions into n orbitals. 

What does it exactly mean, to put parafermions into orbitals? Inspired by the parafermion mode counting derived 
in we present the following rules. While the rules are relatively complicated, they are the only way known to us 
which allows to obtain j "} and, at the same time, to break the configurations into angular momentum multiplets. 

We replace the boxes by the 'positions'. As we move from right to left, each next position carries 1/3 more of the 
z-component of the angular momentum, as in Fig. ^ ranging from — (3ri — 2)/6 to (3n — 2)/6. Thus, we have 3n — 1 
positions. Fig. depicts the positions at n = 2. 
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FIG. 2. Positions where the parafermions are to be put for n — 2 

There are two types of the Z3-parafermions we can put into these positions. We depict them by either the light 
shaded or dark shaded circles. The following rules should be used to put the parafermions into the positions: 

1. The light shaded circle carries L z equal to the L z of its position. 

2. The dark shaded circle carries L z equal to twice the L z of its position. In addition, the dark shaded circle is 
counted as two parafermions. 

3. The dark shaded parafermion can occupy any position, while the light shaded parafermion is not allowed to occupy 
the rightmost or the leftmost position. The number of empty positions to the left of the leftmost parafermion has to 
be 31 if it is dark shaded, or 31 + 1 if it is light shaded, where I is any nonnegative integer, including 0. For example, 
if the leftmost parafermion is the light shaded one, it can occupy the position number 2, or 5, or in general 2 + 31 
counting from the left, and if it is a dark shaded one, it can occupy the position number 1, or 4, . . .. 

4. The adjacent light shaded parafermions are allowed to have 31 empty positions between them, while the adjacent 
dark shaded parafermions have to have 31 + 1 empty positions between them. If a light shaded parafermions is adjacent 
to a dark shaded parafermion, they have to be separated by 31 + 2 empty positions, where I is any nonnegative integer 
including 0. 

Examples of the allowed configurations are shown on Fig. H. This figure depicts the 3 ways one can put 3 



parafermions into 5 positions at n = 2, consistent with (|9|) and (12). Moreover, we immediately observe that these 
configurations form an L = 1 multiplet. 
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FIG. 3. Three possible ways one can put 3 Z3-parafermions into 2 boxes, in accordance with { g } =3. 

It is not very convenient to draw positions all the time, so we introduce the following compact notations. Each 
light circle will be represented by ip s with s equal to L z of the circle. Each dark circle will be represented by <p s . For 
example, the three configurations of the Fig. || can be represented as 

V'iV'oV'^ij ^±0^2; (13) 

(we remember that in accordance with the rule 2, <f> is counted as two parafermions). 
There is only one way to put 6 parafermions into 5 positions at n = 2, 



>1 <P0<, 



(14) 



consistent with { g } = 1 . 

Let us demonstrate the power of the technique by breaking the n = 3 configurations into multiplets. There should 
be 8 positions, carrying L z from —7/6 to 7/6, that is, —7/6, —5/6, —3/6, —1/6, 1/6, 3/6, 5/6, and 7/6. 
It is obvious there is only 1 way to put zero parafermions into these positions. 
According to ( fTo| ) and ( |l2| ) {3} = 10. Indeed, we can count the configurations directly 



L 
L 



«(^i)=5' L »(^*)=|' M^0 = 4' 
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We observe 10 different configurations. Moreover, we observe that L z = ±5/2 occurs only once each, and L z = ±3/2 
or L z = ±1/2 occurs twice each. We immediately recognize the sum of L = 5/2 and L = 3/2 angular momentum 
multiplets. 

Quite analogously, for 6 parafermions at n — 3 we find 



L z (^ps^sTpiTp-iTp-alp-s^ = 0, (16) 

L z (VfV^i^-) = -1, L z (VjVi^I^-) = 1, 

L z (^i^-) = 2, L z (VsV>|^<^) = -2, L,(^^^)=0, 

(02030^1 ] — 3, L z ( 4>zd>3(t> -7 1=1, L z (</>ic/>^30^7 ) = —1, [ d>±4>^3_d>^i ) = —3. 

\ 8 6 6 / \ 6 6 8 / \ 6 6 6 / \ S 6 6 / 

The total multiplicity is 10, in agreement with (|To|) . Counting i z we observe that = ±3 and L z — ±2 occurs in 
one configuration each, while L z — ±1 and L z = occurs twice each. We recognize the direct sum of L = 3 and 
L = 1 multiplets. 

One can further check that it is impossible to fit 9 or more parafermions into these positions. 

The Tables ||, and III in the appendix below contain the numerical data available for the Hamiltonian (Ji|) at 
k = 3 and at the number of electrons fixed at N — 6, N = 9 and N = 12. The data is in the form L versus n. For each 
excess flux quanta n the intersection of nth column and Lth row contains the multiplicity of the angular momentum 
L at given n. 

To reproduce this data using the parafermion counting rules, we need to combine the angular momenta of the 



parafermions with those of bosons, by analogy with (12). This is now trivial to do, since parafermions and bosons are 



distinguishable particles, and their angular momenta combine in accordance with standard rules. 
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For example, at N = 6 we can obtain the multiplets in the following way. First we need to put 9 bosons into 3 
orbitals as represented by the term of ( |l2| ) with F = 0. Their angular momenta are calculated in exactly the same 
way as we calculated the angular momenta for 3 bosons in 3 orbitals in the text directly preceding Fig. |l|. 

Then we have to put 9 bosons into 2 orbitals and combine their angular momenta with the L = 5/2 and L = 3/2 
angular momenta of the parafermions, as represented by the term of ([l2]) with F = 3. 

And finally, we put 9 bosons into 1 orbital. There is only one such state and its angular momentum is 0. Therefore 
the total angular momentum is just L = 3 and L = 1, coming from the parafermion contribution. 

We are not going to go through explicit counting since it is rather standard. The only nontrivial step was the 
parafermion angular momentum contribution, and that was worked out above. We have checked, however, that the 
result reproduces the numerical data given in the Table [jj. In fact, we have verified all the data in all three Tables, 
expect the results for n = 5 at N = 6. In particular, we found that the parafermion multiplets at n = 4 are the 
following. At F = 3, there is one each of L = 0, 2, 3, 4 states giving the multiplicity of 22 (compare with (pp|)). At 
F = 6, there are two each of L = 0, 2,4 states, and one each of L = 3,6 states giving the overall multiplicity of 50. 
For F = 9 it is one of each L = 0, 2, 4, with a total multiplicity of 15. And for F — 12 there is one L = state. 

It would be nice if we could obtain the parafermion multiplicities { F } in a more closed form rather than by having 
to count configurations all the time. The most closed form that we are aware of can be obtained with the help of the 
equation for the partial partition functions derived in fl5| . 

Wc define the polynomials Yi(x) in the following way. 

Y l+1 (x) = xY l+l (x) +x 2 Y l+1 _ (x) + (1 - x 3 )^), (17) 
with the initial conditions Y_i = 0, Yq = 1, and Yi = 0. Then we claim that the following expansion is valid 

Y »w= £ ( 18 ) 

fc=0,l,... 

It is in fact possible to derive ( |l7| ) directly from the counting rules presented above, see fl5f| . 
It is not hard to check that 

Yi = 1, Y 2 = 1 + 3.t 3 + x 6 , Y 3 = 1 + 10.x 3 + 10x 6 (19) 



in agreement with (J9j) and (10). We can continue further, 

y 4 = l + 22.T 3 + 50.x 6 + 15x 9 + a; 12 , Y 5 = 1 + 40.x 3 + 168a; 6 + 140x 9 + 28x 12 . (20) 

These multiplicities are also in agreement with the available numerical data for the degeneracies at various flux quanta 
n. 

We note that according to the sum of the multiplicities at fixed n should give us Fibonacci numbers F 3n ^2, 



£ { 3 ;M„-, (2D 



.3k 

fc=0,l,2,... 



That is indeed true, since this sum is equal to Y n (l), and by substituting x = 1 into ( [l7| ) we obtain 

Y l+1 (l)=Y l+i (l)+Y l+i (l), (22) 
which is the defining relation for the Fibonacci numbers. 



III. CONFORMAL FIELD THEORY AND THE COUNTING RULES 

The rules we have just presented were in fact deduced from the parafermion conformal field theory. Here we would 
like to explain how conformal field theory should be employed to derive these rules. 

While we were interested in parafermion quantum mechanics, conformal field theory, as a field theory, gives us 
the counting rules in the second quantized formalism. It is not so hard to read statistics off the second quantized 
formalism. 

Let us concentrate on the Majorana fermions as an example. We know, of course, that there are (^) ways to put 
F fermions into n boxes, in accordance with Pauli exclusion principle. This can be obtained also from the conformal 
field theory of the Majorana fermions. 
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Expanding the Majorana fermion %j) n in terms of modes we get 

#s) = X>«*- n -*, (23) 

n 

with n going over half integers 1/2 + Z. The modes ?/>„ obey the anticommutation relations, 

^n^m + 1pmi>n = <Wm,0- (24) 

All these are well known facts of conformal field theory. 

Since the Majorana fermion defined at every point in space and time cannot be infinite when acting on the vacuum, 
ijj(z)\0) for any z including z = 0, it follows that ip n \0) = for all n — 1/2, 3/2, . . .. We say that the modes of 'tp(z) 
with positive index act as annihilation operators. The modes ij} n with n = —1/2, —3/2, . . . act as creation operators. 
We can use them to create new states. A generic state will look like 

iji- nF ip- np _ 1 ...ip^ ni \0), (25) 



with all the rij being negative half integers. Note that some of the states in (25) are not linearly independent and 
can actually be transformed into each other by repeatedly using the anticommutation relations (p4|). We can select a 
linear independent subset of ( pq ) by ordering all the n, say making them increase from left to right, and making sure 
neither of them are equal to each other, thereby making the state zero according to ^„ = 0. 

Now observe that there is the following correspondence between the states (Eq) and the angular momentum wave 



functions. Consider the states (Eq) together with the relevant powers of the coordinates, as in 



Ip-nplp-nF-! ■ ■ •V'-nJO) z" F 2 z'% " 2 . . . Zp 2 . (26) 

Let us now sum (^) over all the states accessible via the anticommutation relations (0). In other words, we want to 
sum over all the permutations a of the numbers rix, n 2 , • ■ ■, np. Obviously we obtain 

W-<T(n F )V-<T{n P -i) •••V-o-(ni)IO) z l z 2 2 ...z p 2 = (27) 

ip- nF ,il>- nF _ 1 ...ip- ni \0) ^sign(a) z{ { 



^cr(n F )-i cr(njr_i)-5 cr(ni) — 5 

z 2 . . . ;Zp 



One immediately recognizes the totally antisymmetric polynomials of the sort discussed in |14|. The angular momen- 
tum L z computed on these polynomials can be defined, up to an additive constant, as generated by the operator 
to give rii — F/2. These totally antisymmetric polynomials are in fact the wave functions of the fermions in the 
first quantized formalism! 

If we consider all possible states (p7|), with the restrictions tin < n, n being some integer, there is going to be 
exactly (^) of them which are the multiplicities in ([ll]). Now we see how these states break naturally into angular 
momentum multiplets. One can check that the multiplets we obtain in this way matches the numerical data for the 
Pfaffian state Q. 

What we just did so far looks like a complicated and not very intuitive way of rederiving the results of fll4| . However, 
as we move on to the parafermion states, conformal field theory becomes the only consistent and dependable way to 
construct the angular momenta multiplets. 

k = 3 parafermion conformal field theory contains two fields of dimension 2/3, ip^(z) and ip^(z) |T^ ]. Expanding 
them in terms of modes we obtain 

^) = Ev^— f , v (2) (*) = £^ 2) ^ s - (28) 

n n 

It is obvious that the modes ipn^ and ip^ have to annihilate the vacuum if n > —2/3. Applying the modes with 
ri < —2/3 to the vacuum we can in fact create parafcrmionic states. Not all the states we create in this way are 
linearly independent. To check which ones are independent we have to employ the generalized commutation relations 
of the sort derived in These relations are very complicated. Fortunately for k — 3 the set of independent states 
have already been derived in fL5| |. Here for completeness we are going to quote the answer. 

The full set of linearly independent states is created by applying the modes of the parafermions of the first kind 
<f>s = ipg or a certain combination of the modes 4>^ = ips ip _2- We do not need to use the modes of the 
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parafermion of the second kind as its modes create states which are linearly dependent on the states created by 

4 1] and 4 2) - 

The allowed states have the form 



with the spacing specified as 



ej-.-^^IO), (29) 



if ij+i = 1, ii = 1 then s;+i - si =m+\, (30) 



if ii+i = 2, k = 1 then si +1 - si = m 
if = 1, ii = 2 then s; + i — si = m 
if = 2, = 2 then s i+ i - si=m 



3' 
2 

3' 
4 

3' 
2 

3' 



where m is any nonnegative integer number including 0. 

We note that these states are in one to one correspondence with the counting rules of the previous section. In fact, 
this is how the counting rules should be derived. 

We can go slightly further and conjecture a way to derive the parafermion wave function. We need to consider the 
sum over all the states which can be obtained from one of the states (E9h by the generalized commutation relations 



£ ^...^f)|0) z-"...^ 1 ^, (31) 

dependent states 

with different j being either 1 or 2. By using those generalized commutation relations we can in principle bring the 
sum to the form (compare with (|27|)) 

^...^i^JO) f(zi,...,z N ). (32) 

We interpret f(zi, ■ ■ . , z/v) to be the wave function of the parafermions! 

While finding / explicitly is beyond the scope of this paper, we note that it indeed combines the single particle 
states of the form z s into the multiparticle wave function obeying the right particle exchange properties. 

One should also in principle be able to derive / from the quasihole wave functions in direct analogy with the 
corresponding derivation for the Pfaffian state |13|Jl4| ]. This also lies beyond the scope of this paper. 

At the end we would like to note that for k > 3 the mode counting rules can also be derived by using the generalized 
commutation relations of p2| . Some partial results at k ~ 4 and k = 5 are presented in pl| . 



IV. THE CORRELATION FUNCTIONS OF THE PARAFERMION CFTS 

All the preceding discussion of this paper was based on the assumption that the wave function (^) can be generated 
as a correlation function of the parafermions, in the sense of j|. 

While it was conjectured in that (^) is such a correlation function, no proof was found. It is important to prove 
this relationship, otherwise, our manipulations with parafermions lose their relevance. 

In this section we would like to present a proof that the wave function found in Q], which is the zero energy 
state of the k + 1 particle interaction Hamiltonian (|]) , is indeed the correlation function of a parafermion conformal 
field theory. For this purpose we recall that a parafermion conformal field theory consists of k — 1 fields iph 
I = 1, . . . , k — 1 with the operator product expansion 

M^Az') KQ^z-zr 2 ' 7 '^ + p liV (z - z')dxPi +v + ...), (33) 

where c and j3 are certain numbers - structure constants. The index I of the fields tpi in (|3^) has to be understood 
in the mod k sense, ipk+i = V'z- Additionally, ?/>o is identified with the unit operator. In that case, dtpa vanishes and 
therefore, when I + I' = mod fc, the linear term in ([33]) proportional to j3 vanishes. 
According to the conjecture of M the correlation function of N fields ipi is equal to 
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(ipi(z 1 )...ip N (z N )) = — ^ (34) 

\Ai<M - z i) k 

with $ given by (||). 

To show that the right hand side of is indeed consistent with (33), we will glue k — 1 parafermions together 
to obtain the correlation function with tpk-i- Then we will glue ipk-i and ipi together and show that the correlation 
function is consistent with the fact that 

^-i(z)ViOz') <x(z- z 1 )- 2 ^^ + ... (35) 

with a vanishing first derivative dipo. 

The proportionality sign in ( |3^ ) and throughout this section means that an unimportant numerical constant has 
been dropped. 

We start with taking the limit z\ — > z 2 . In that limit the right hand side of (|34|) is proportional to {z\ — z 2 )~ 2 ^ k ■ 
This is indeed consistent with the operator product expansions 

M Z )M Z ') <*{z- zTh2(z') + ... (36) 
By multiplying ( |34"| ) by (z\ — z 2 ) 2 / k and taking z\ — > Z2 we arrive at the following correlation function 

I L t \ i. ( \ it \\ ^(z 2 ,z 2 ,z 3 ,...,z N ) 

{lp2{Z2m{Z3)...1pl{z N )) oc — -j— : -5-. (37) 

n i > 2 ( ;2 2 - ZiP [[ 2<1< j<n(^- z i) k 



As Z2 approaches z 3 the right hand side of (|3^) has just the right singularity (z 2 ~ z 3 )~ A / k , matching the singularity 
of ([53]) as V^2 approaches "01 ■ Therefore, we continue this process further and take z 2 — > z 3 . By doing so we indeed 
recover the correlation function with the field tp 3 , 

I I ( \l ( \ ll \\ 1 f{z 37 z 3 ,z 3 ,z 4l ...,z N ) 

(il) 3 {z 3 )ip 1 {z i )...%l) 1 {z N )) oc — - T — -T- (38) 

lli>3(^3 - Zi)k [[ 3<i<j < N {Zi ~ Zj)h 

It is clear at this point that as we continue 'gluing' parafermion fields together at some point we will arrive at 

/ , / \l f \ I ( \\ ^(Zk-l, Zk-l, ■ ■ ■ , Z k -1, Z k ,Z k+ l, . . . , Z N ) 

{ipk-i{zk-i)Vi{z k ) ■ ■ ■ipi{z N )) oc ^rry 7- (39) 

n^fc-iOfe-i-^) k u k -i<i<j<N( z i - z j) k 

At this stage we should be careful. By taking Zk-i to z k we should be able to recover the identity operator ?p k = iJjq. 
Let us check that this is indeed the case. The main singularity of ( p9| ) as — > z k is (z k _i — z k ) ~ which indeed 
matches (|3^) and therefore, 

III \ I ( \ II \\ ^(z k ,z k ,...,z kl z k+ i,z k+2 ,...,z N ) 

{tpo{zkm{zk+i) ...tpi{z N )) oc — — — -5-. (40) 

Ui>k{ z k ~ z i) 2 Uk<i<j<N( z i - z 3) k 

By its definition, the identity operator does not depend on its position, and therefore the right hand side of ( |40| ) 
should not depend on Zk- To see that, let us recall that in [Q the following theorem was proved. The wave function 
^ for N particles whose first k particles live at the same point Zk can be expressed in terms of the wave function for 
N — k particles in the following way, 

W(zk, Zk, ■ ■ ■ , Zk, Zk+i, ■ ■ ■ , z N ) oc ] [ (zk - z l ) 2 *(zfe+i,Zfe +2 , • • ■ ,z N ). (41) 

k<i<N 

Substituting © to © we see that 

1 1 1 \i 1 \ 11 \\ *( z fe+i' z k+2 , ■ ■ ■ , zn) ,. nS 
(ip {Zkm{zk+i) ■ ■ >ipi{zN)) oc — -f, (42) 

llk<i<j<N\ z i ~ z j) k 

that is, is indeed an identity operator. The correlation function is insensitive to its insertion. That completes the 
first part of our proof. 
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If we continue to expand (39) in powers of Zk-i — Zk we observe that we should reproduce the operator product 
expansion ([53]) with I + 1' = mod k. In particular, we must see that the linear term of that expansion vanishes, 
flk-i.i = 0. If this term didn't vanish, that would mean that in addition to the parafermion fields, the conformal field 
theory whose correlation function is given by ( p4| ) has a dimension 1 operator J which generates a U(l) symmetry 
(see [Q). Such an operator should be absent in the parafermion theory. Let us check that it is indeed absent. 

Continuing the expansion of (p3) in powers of Zk-i — Zk we find that the linear term is obviously proportional to 



d 



dz k - 



^(z k -l,Z k -l, ■ 



. Zk-l, z k, z k+l, ■ 



IL>kOk-i 



2(fc-l) 



,Z N ) 



(43) 



To show that this term is zero it is sufficient to show that 



-^(Zk-l, Zk-l, ■ ■ ■ , Zk-l, Zk, Zk+l, ■ ■ ■ , z N ) 



dz k - 



2(fc — 1) 



\i>k 



Zk - Zi 



^(z k ,z k , ■ ■ .,z k ,z k +i, ■ ■ - ,z N ). 



(44) 



We can further simplify the equality ( |44j) by noting that since ^ is a symmetric function of its arguments, it is enough 
to differentiate it with respect to just one variable, 
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- — ^ (Zk-l, Zk-l, 

dz k 



. Zk-l, Zk, Zk+l, ■ 



,z N , 




^{z k , z k ,..., z k , z k+ i, z N )- (45) 



To prove ( |45"| ) we again recall the definition of ^(zi, . . . , zn) as given in (^). To construct it, we break all its coordinates 
into N/k clusters of k coordinates and then sum a certain expression, written in terms of these clusters, over all the 
permutations of N particles. Let us first look at the right hand side of (ft5|). In accordance with the theorem proved in 
0, the only terms which contribute to ^(zk, ■ ■ ■ , Zk, Zk+i, ■ ■ ■ ,zisr) in the sum over permutations (||) are those where 
all the Zk belong to the same cluster. This is how (|4l]) could be derived. Let us now make an assumption, which we 
will justify later, that the only terms which contribute to the left hand side of (ff5|) when we substitute (H) for \& are 
also those where all the Zk-i and Zk belong to the same cluster. 

It is possible to convince oneself that after some algebra the sum of these terms can be reduced to 



4> 



ee n 

r=l P Q<j<N/k 



(z k — Zp( r+ kj))(zk — Zp( r+ kj+s(r))) 
, [Zk-l — Zp(r + kj)){Zk-l — Zp( r +kj+s(r))) , 



Yl (zk-i-zif Yl Xj+i,i+i, 



(46) 



k<i<N 



0<j<KN/k 



where P(n) gives a permutation of the integer numbers k + 1, k + 2, . . ., N. And s(n) = n + 1 for < n < k with 
s(k) = 1. x are t ne expressions (||). Note that \ do not depend on Zk-i and Zk- Applying the derivatives as in ( |45| ) 
we get 



d 



dz k 

k 

EE E 

r=l P 0<j<N/k 



(47) 



1 



1 



1 



,k<i<N 



Zk Zi 



Zk — Zp( r+k j) Zk Zp( r+k j+s(r)) 



4>l 



Yl (z k - Zif Yl xj+i,i+i = 

0<j<KN/k 



k<i<N 



The last line in ( ^7| ) follows from the fact that the summation over r completely symmetrizes the sum in the second 
line of fl47| ) over all the permutations P. 

Therefore we have proved ( |4^ ) on the assumptions that the only terms that contribute are those where the first 
k coordinates of '5 belong to the same cluster. To see that it is indeed so, let us recall again that according to the 
theorem proved in all such term should vanish as Zk-i approach Zk- Since they are polynomials they vanish at 
least as Zk — Zk—i, or perhaps even faster. If they vanish faster, their derivative with respect to Zk with the setting 
Zk-i = Zk after differentiation is definitely zero. If, on the other hand, they vanish linearly, then we can argue that 
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for each term vanishing as Zk — Zk-i there is another permutation of the coordinates different from the first one 
by exchanging exactly the two coordinates in that difference. This term will be proportional to z^-i — Zk- After 
differentiating and setting Zk — Zk-u these terms will cancel each other. 

This concludes the proof that the wave function conjectured in to be the correlation function of the parafermions 
is indeed the correlation function of the parafermions. This in fact allows us to use the parafermionic conformal field 
theory to describe the ground states of the hamiltonian (0), as we did throughout this paper. 
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APPENDIX A: NUMERICAL DATA 

In this appendix we present the numerical data for the degeneracy of the angular momenta multiplets at various 
excess flux quanta n. Columns of the tables are labeled by n. Rows are labeled by the angular momentum L. The 
intersection of the nth column and Lth row gives us the degeneracy of the angular momentum L multiplets at a given 
n. The three tables represent the angular momentum multiplets at the total particle number N — 6, N = 9, and 
n = 12. All the data is at k = 3. 
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TABLE I. N = 6. 
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TABLE II. N = 9. 
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TABLE III. TV = 12. 
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